Integer Quantum Hall Effect in Double-Layer Systems 
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We consider the localization of independent electron orbitals in double-layer two-dimensional elec- 
tron systems in the strong magnetic field limit. Our study is based on numerical Thouless number 
calculations for realistic microscopic models and on transfer matrix calculations for phenomenolog- 
ical network models. The microscopic calculations indicate a crossover regime for weak interlayer 
tunneling in which the correlation length exponent appears to increase. Comparison of network 
model calculations with microscopic calculations casts doubt on their generic applicability. 
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I. INTRODUCTION 

The integer quantum Hall effect is generally well un- 
derstood in single-layer two-dimensional electron systems 
(2DES's) which are sufficiently disordered that interac- 
tions do not play an essential role and are in a field suf- 
ficiently strong that Landau level mixing does not play 
an essential role. In this limit, single-electron orbitals 
are localized except at a critical energy E c near the cen- 
ter of each disorder broadened Landau level. For Fermi 
energy E F = E c , theoryLTQ predicts that {cr% x ,cr% y ) = 
(1, 2n + l)e 2 /2h whereas on the Hall plateaus (Ep ^ E c ) 



[p xx 



(0, n)e /h where n is the number of ex- 



tended state energies below the Fermi level. As the criti- 
cal energy is approached, the localization length for elec- 
trons at the Fermi level is expected to have a power- 
law divergence, £ ~ \Ep — E c \~ v , and v, the correla- 
tion length exponent, is expected to be independent of 
n. It is believed that-ihe transition is well described by 
quantum percolatioroj models and semiclassical calcula- 
tionsQ have estimated the correlation length exponent to 
v = 7/3. This picture has been corroborated by a large 
number of thorough numerical studiesat3 which are in 
agreement with theoretical predictions for a xx and a xy . 
Localization properties and the divergence o£-th|e-|local- 
ization length have been studied extensively ll3~E3 with 
the most recent estimate oi-jthe correlation length expo- 
nent being v = 2.35 ± 0.03B3. On the experimental side, 
measurements of the width, AB, of the peak of p xx as 
well as (dpxy/dB) -1 , both predicted to scale with tem- 
perature as T l / zu , yield values 1/v = 0.42 ± 0.0400 
with z assigned to be 1. Higher derivatives of p xy yields 
exponents^ of nv in agreement with scaling theories of 
the transition between Hall plateaus. (Experiments in. 
the fractional quantum Hall regime find similar valuesEj 
for this exponent.) Recently thf? dynamical critical expo- 
nent, z, has been measurecSc3 to be z = 1. 

In this paper we report on a numerical study of the lo- 
calization properties of single-electron orbitals in double- 
layer two-dimensional electron systems. This work is mo- 
tivated by recent experiments hinting at changes in local- 



ization propertieSr-wheiar-two different Landau levels are 
nearly degenerateBHrE3'EZrE£l, by growing interest in the 
conditions necessary for the occurrence of the quantum 
Hall effect in three-dimensional electron systems, and by 
the need for improved understanding of the disappear- 
ance of the quantum Hall effect at weak magnetic fields 
in high mobility samples. In each case, we believe that 
double-layer quantum Hall systems offer advantages for 
both theoretical studies and for the experimental studies 
which wc hope to motivate. 

In single-layer two-dimensional electron systems, local- 
ization properties appear experimentally to be changed 
when the exchange enhanced spin-splitting between. Lan- 
dau levels with the same orbital index collapses. EJ The 
interpretation of these experiments is confused by un- 
certainties involved in modeling the spin-orbit disorder 
scattering necessary for mixing the two Landau levels and 
by the the apparent importance of interaction effects in 
controlling the degree of mixing. The interaction com- 
plications are not so troublesome in double-layer systems 
and, in addition, the degree of mixing between Landau 
levels in separate quantum wells can be controlled by ad- 
justing the strength of the barrier separating the wells 
or by adding an external bias potential which moves the 
double- layer system off balance. 

In high-mobility two-dimensional electron systems 
(Eft /ft, >> 1), the quantum Hall effect appears to be- 
come unobservable in practice once Landau level mixing 
by disorder becomes strong, i.e. once uj c t is of order one. 
(Here lu c = eB/m*c is the cyclotron frequency.) The loss 
of an observable quantum Hall effect in these systems 
appears to be associated with a dramatic increase in the 
localization length in the middle of the Hall plateaus, 
rathpjthan with the 'floatation' of extended state ener- 
gies EIlEJ which occurs in more strongly disordered sys- 
tems. It seems likely that the same dramatic increase in 
localization lengths on Hall plateaus will occur in double- 
layer systems when the Landau levels in the two-layers 
are strongly mixed. The ability to systematically control 
the number of Landau levels which are mixed motivates 
working with double and multi-layer systems. 

Since much of the physical picture underlying the 
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quantum Hall effect is specific to two dimensions, it 
was not initially clear that it was even possible to ob- 
serve quantized plateaus in three-dimensional systems. 
Early experimental work, focused on widely separated 
2D layersEa, found the IQHE with quantized resistiv- 
ity Pxy = h/N z ie 2 where N z is the number of quantum 
wells, consistent with parallel conduction in many quan- 
tum wells. Stormer et alS3 performed experiments on 
coupled GaAs superlattices with 30 periods, where the 
dispersion relations explicitly showed three-dimensional 
effects in zero magnetic field. Measurements of the re- 
sistivity tensor showed that p xy again was quantized as 
h/nie 2 despite the 3D nature of the system, but it was 
found that n ^ N z . This was later explained in terms of 
band bending which raised the energy of states in quan- 
tum .wells close to the sample surface above the Fermi 
level.c3. Subsequent experiments on a 200-period super- 
lattice seem to confirm this picture since only a_Bxed 
number of apparently empty quantum wells occuilJ. In 
later work Stormer et a/M3 demonstrated that o~ zz also 
has deep minima on quantum Hall plateaus. Significantly 
for the physics addressed here, no Hall plateaus with in- 
termediate integer indices were observed. 

The phase diagram of disordered three-dimensi«riaJ| 
systems in a magnetic field has been investigatedE3~c£l 
and the possibility of a metal-insulator (MI) transition 
has been pursued. The study of the quantum Hall ef- 
fect in three-dimensional systems may also be relevant 
to quasi-one-dimensional systems such as the Bechgard 
salts which ioiun a spin density wave (SDW) state in a 
magnetic fieldE! In such a state ttuiJiall conductance is 
also quantized- Recent theoreticaOO as well as exper- 
imental workEJ supports a picture in which a complex 
phase diagram arises from the interplay between SDW 
formation and the quantum Hall effect. 

Theoretical model calculations have examined the An- 
derson transition in a strong magnetic field for thixp* 
dimensional systems. Transfer matrix aalculationsSEHI 
and recursive Green's functions studiesEj have shown 
that the localization length at the mobility edge diverges 
with an exponent v = 1.35 ± 0.15E3 similar to what is 
found in the absence of a magnetic fielcO. Calculations 
performed on network modelsEj indicate an exponent of 
v = 1.45 ±0.25 in, perhaps, surprisingly good agreement. 
Note that the exponent v is significantly smaller than the 
approximately 7/3 found at the ID quantum Hall tran- 
sitions. 

In our studies of double-layer systems we assume that 
the disorder potentials in the two layers (see Fig. [I]) 
are uncorrelated. As sketched in Fig. [l] it is crucial to 
make the distinction between this form of disorder, which 
we shall refer to as uncorrelated disorder, and the form 
where the disorder is identical in the layers (correlated 
disorder) since the latter has a much smaller effect on 
the localization length. Related work on the double- 
layer system has previously been done by Ohtsuki et 
al£3. Our study has clear analogies to previous work 
on the spin-degenerate (where the two spin levels are not 



resolved) and spin- resolved transitions. Experimentally, 
the spin-degenerate transition has beep^investigated in 
several different studies. Wei et. al.c^l observed that 
AB as well as (dpxy/dB)^ 1 both behave as a power-law 
in T with a much smaller exponent than obsep/ed for 
the spin-resolved transitions. Later experimentsEZI found 
AB ~ T 0,21 . Experiments on the frequency, (/), depen- 
dent conductivity have shown that a xx peaks broadens 
as {AB) ~ p with 7 = 0.41 ± 0.04 for the spin split 
transition and 7 = 0.20 ± 0.05 for the spin-degenerate 
caseE3. This, tantalizing effect was explained in terms of 
an unstable critical point so that the enhancement of the 
exponent-, is understood as an artifact of crossover phe- 
nomenaEZI (See Fig. |J). However, other workEJ has pro- 
posed the possibility of an effectively stable fixed point 
due to a disorder induced destruction of the exchange en- 
hancement of the electron g-factor. In the double-layer 
system this would correspond to a finite critical coupling 
t c needed to see the symmetric and antisymmetric state 
appear. Network model calculationsElfEa find that the 
universality class is unchanged in the presence of strong 
Landau level mixing between the polarized Landau sub- 
bands. However, it is presently not clear to what extent 
these models actually describe Landau level mixingti^Ei 
Thej-efjEect of spin-orbit scattering has also been consid- 
erecErEll again leading to crossover phenomena without 
a change of the universality class. In the language of the 
double-layer systems these studies consider completely 
random tunneling between the two layers and in most 
cases the disorder is strongly correlated between the two 
layers. Since the experimental situation for double-layer 
systems is much closer to a weak and uniform tunneling 
between the two layers, which was not considered pre- 
viously, we investigate this limit in detail. In the case 
of the spin-degenerate spin-resolved transition the con- 
trolling parameter is the electron g-factor which is not 
readily tuned experimentally. However, for double-layer 
systems, it should be experimentally possible to tune the 
coupling between the two layers by tilting the magnetic 
field. This is a very convenient circumstance since it is 
(conceivably) possible to investigate if a finite coupling 
t c is needed in order to split the Landau level subbands 
corresponding to the symmetric and antisymmetric state. 
Assuming that the interesting interaction effects associ- 
ated with tilting the magnetic fieldEiliEl can be ignored 
in the more disordered samples to which the considera- 
tions of the present paper would apply, the only effect of 
tilting the magnetiCriield in a sample with uncorrelated 
disorder is to reduceEil the effective tunneling parameter: 

t* = icxp(-(dtan(6>)/2f) 2 ) (1.1) 

where d is the layer separation, 2n£ 2 B± = $0 with B± 
the component of the magnetic field perpendicular to the 
sample, and 8 the angle by which the magnetic field is 
tilted away from the normal to the layers. 

Our paper is organized as follows. In Section || we 
discuss several different arguments, all at the semiclas- 
sical level, which indicate that even very weak coupling 
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between the layers can give rise to dramatic increases in 
localization lengths and possibly change the apparent lo- 
calization length exponent. In Section III we present the 
results of a numerical study for a realistic microscopic 
model of a double-layer system. Localization properties 
are discussed in terms of the Thouless number g s . Sec- 
tion IV is concerned with network model suitably defined 



to describe the weakly coupled double-layer system. Us- 
ing transfer matrix techniques we calculate the reduced 
correlation length. Finally in Section ^ we present our 
conclusions. 



II. SEMICLASSICAL ARGUMENTS 

When the disorder potentials in both layers are smooth 
on microscopic length scales electronic orbitals are lo- 
calized along equipotentials. The quantum percolation 
theory of the integer quantum Hall effect is based on a 
theory of percolating equipotentials supplemented by the 
possibility of quantum tunneling between equipotentials 
near saddle points of the disorder potential. Semiclassi- 
cally, the 'E cross B' drift velocity of an electron along 
an equipotential is proportional to the local electric field. 
For a Landau level of width T the typical electric field 
is T/eA where A is the correlation length of the disorder 
potential. Accordingly the typical drift velocity is 



V dr = 



cT 

zWR' 



(2.1) 



We are interested in the influence of tunneling between 
the layers on localization properties. Tunneling intro- 
duces a new typical length scale into the physics of the 
system, the drift length 
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(2.2) 



Idr is the typical distance an electron drifts along an 
equipotential between tunneling events. Typical equipo- 
tentials are closed paths with a perimeter which we 
can for present purposes associate with the localization 
length £ whichrdiverges at a critical energy within each 
Landau levelSO. For weak tunneling is long. When 
Idr is larger than the sample size, tunneling will have no 
effect on even the most extended orbitals in the system 
and we should expect negligible changes due to tunneling 
in our finite size system calculations. When ld r is smaller 
than £, electrons will typically tunnel before completing 
a closed orbit around an equipotential. Note that this 
condition is satisfied at smaller and smaller t values as 
the center of the Landau level is approached. When an 
electron tunnels to the other layer, it will move along an 
equipotential of a statistically independent smooth disor- 
der potential. When it later returns to the original layer, 
it will not in general return to the same equipotential 
contour from which it departed even neglecting tunneling 
events near saddle points within a layer. Related ideas 



have been discussed in Ref. [30| In our view the pos- 
sibility that localization physics is qualitatively altered 
by tunneling between the layers deserves serious atten- 
tion. The situation becomes simple again only when ld r 
is smaller than A so that an electron tunnels many times 
before its local potential profile changes. In this limit, 
electronic eigenstates will be symmetric and antisymmet- 
ric combinations of the individual layer eigenstates and 
the effective disorder potential will be the mean of the 
independent disorder potentials in the two layers. Note 
that the limit of infinite system size and the limit of van- 
ishing tunneling amplitudes are not interchangeable. In 
the thermodynamic limit, tunneling will always be im- 
portant for those states near the critical energy which 
have a localization length larger than ld r - 

In Fig. ^ we show the flow diagram proposed in Ref. |27]. 
The controlling parameter that describes the flow from 
the unstable fixed point, P, to the more conventional pic- 
ture (solid lines) is the electron g-factor or the tunneling 
parameter, t, for the double layer systems. Note that the 
unstable fixed point is at + since two uncoupled layers 
(t = 0) cannot lead to new critical behavior. In Fig. ^ we 
have left out the flow around the unstable critical point 
P leaving open the possibility of a new unstable critical 
point at a finite t c . Such an unstable critical point would 
imply that the transition from the v = phase would be 
directly into v = 2 phase without an intervening v = 1 
state. This scenario has recently been discussed by Tikof- 
sky and KivelsonEJ and recent experiments in the strong 
disorder weak field limit could Jae interpreted as lending 
support to this possibilityESEZTEJ. 

If, as suggested in Ref. [27], the enhanced exponent seen 
in the experiments on the spin-degenerate transition is 
indeed due to the presence of an unstable fixed point one 
might ask whv— the crossover exponent should be 14/3. 
Polyakov et a/J!3 have proposed a crossover form for the 
correlation length based on the assumption that the ef- 
fective correlation length is the square of the correlation 
length in the absence of any coupling. We now give a 
brief argument, somewhat speculative and heuristic, in 
the spirit of the semiclassical calculation of Ref. for 
why this could be the case. In the absence of any tunnel- 
ing to the other layer we will at a given energy E, relative 
to the middle of the Landau level, have states localized 
on a equi-potential contour £q(E). We then introduce a 
very weak coupling, t, to the other layer. Following the 
discussion in the beginning of this section we expect that 
for £o(E) <C ldr there will not be sufficient time to scat- 
ter into the next layer before the electron self- interacts. 
Thus a non-zero t will not affect sufficiently small £o(E). 
However, when £o(E) ~ ldr, but still in the limit («T 
the electron will scatter a number of times, N s , along the 
final path £(E). Naturally, 



N s ~£ {E)/l dr . 



(2.3) 



Furthermore, in the limit t -C T we expect the electron to 
be scattered among different orbits in different layers at 
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approximately the same energy and spatial extent, £o(E) 
(see top panel in Fig. |l|). We then find that 



\E\ 



(2.4) 



In a given orbit of extent £q(E), which makes up a part 
of the final path of length, £ this argument assumes only 
a few scattering events into different orbits. That is, we 
implicitly exclude events where the electron immediately 
is scattered back into the same orbit. Clearly when t be- 
comes large many scattering events will occur and we ef- 
fectively form symmetric and antisymmetric states. This 
will occur when t > T or equivalcntly 



Idr < A. 



(2.5) 



Thus we expect an enhanced correlation length and an ef- 
fective doubling of the correlation length exponent when- 
ever 



A < l dr < f. 



(2.6) 



III. EXACT DIAGONALIZATION RESULTS 

We now turn to a discussion of our exact diagonaliza- 
tion results for double-layer systems. We shall work ex- 
clusively in the lowest Landau level approximation. Since 
we are considering a finite system of dimensions L x ,L y 
we want to impose periodic boundary conditions. In this 
case one uses the folkwdng set of basis functions for the 
lowest Landau levellirO: 



<f>7(x,y) 
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(3.1) 



Here, Xj jS = j2n£ 2 /L y + sL x and I — ^J%c/2eB is the 
magnetic length and j runs from 1 to N$ = L x L y /2ir£ 2 
where is the number of flux quanta or the degeneracy 
factor of the lowest Landau level. Since we describe each 
quantum well by this set of lowest Landau level wave- 
functions we include the dummy index n to denote the 
different layers. It is easy to see that the individual terms 
in the infinite sum all are invariant under y — > y+L y . The 
sum over s makes the wave- function invariant up to phase 
factors also under the transformation x — + x + L x . We 
model the randomness as 5-function scatterers at random 
position with random sign. 



27rA£ 2 5 n ,„^0r(r p )^'(r p ) 
P =i 

— tSi y j(S n ^ n -i + #n',n+l)) 



(3.3) 



where A contains the random sign of the (5-function scat- 
terers and their stren gth. n p is the number of scatterers. 



The first term in Eq. (3.3) is the potential energy in each 
of the wells while the second describes the tunneling be- 
tween the two wells in the tight-binding approximation, 
with t the tunneling parameter. We shall always use peri- 
odic boundary conditions in the z-direction so that layer 
N z + 1 is identical to the layer 1. We chose |A| so as to 
fix the width of the Landau levels to be of-the order of 
1 in the self-consistent Born approximatiorQ, Tscba, in 
our units we have 



r SC BA = 2|A|y / n p /iV^, 



(3.4) 



and we therefore chose |A| = (l/2)y/N^/n p . With this 
choice we have effectively chosen our energy scale and 
all of our results have energy, E, in units of Tscba- F° r 
different sizes we always keep |A| and therefore thc-j-a- 
tio np/Nff, constant. From exactly solvable modelsO it 
is known that the density of states exhibit peculiarities 
when n p /N<p < 2 we have therefore chosen always to work 
with rip/Ntf, = 3. We have explicitly checked that in the 
case of a single layer we find, with the above mentioned 
definitions that the calculated density of states is well 
described by the following approximate formula 



2ME) = sfl n f e-m\ 



(3.5) 



with, r = 1. In the case where we do not include the 
random sign of the scatterers we have also checked that 
the density of states corresponds to the exact result of 
Ref. [7^. Note that we use the standard definition of p(E) 
which integrates over E to l/(2nl ). 

We note that with our definition of the potential the 
variance of V: 



< V(r)V(r') >= 6irf\ 2 5{r - r 



v 2 £ 2 



S(t- 



(3.6) 



With our choice of A we then have v = y/n/2. We could 
equally well have chosen our energy scale by making the 
choice v = 1, and thereby fixing A. 



V(r) 



2ir\£ 2 S ( r 



(3.2) 



Although this is not a realistic model of real randomness 
it is generally believed that the form of the randomness 
is irrelevant, see however Ref. |l2|. The Hamiltonian can 
then be written 



A. Computational Method 

The Thouless numberEl, g s {E), is defined as the ab- 
solute value of the shift of a given energy level, under 
a change in boundary conditions from periodic to an- 
tiperiodic, lAi^l, multiplied by the total density of states, 
N(E), 
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g a (E) = N(E)AE. 



(3.7) 



Here N(E) integrated over all energies is the total num- 
ber of states, NzNcj,, for an N z — layer system, or in other 
words, N(E) = L 2 N z p(E). Clearly extended states are 
much more sensitive to a change in the boundary con- 
ditions than localized, and g s (E) therefore measures the 
stiffness at a given energy, E. In our calculations we 
change the boundary conditions in the y-direction in the 
individual layers simply by per form ing the transforma- 



tion Xj^ s — > Xj+i/2 >a in Eq. ( |3.3| ). It is known that 
in the absence of a magnetic field the Thouless num- 
ber is related to the longitudinal conductivity, cr xx (-E)Ea. 
For a recent account of Thouless number calculations see 
Ref. |59| and references therein. We follow Ref. ^ in deriv- 
ing a scaling function for the integrated Thouless num- 
ber. If we assume that the correlation length diverges 
as £ ~ we can write a finite-size scaling form for 

9s(E): 



g s = ~g s {EL l ' v ), 



(3.8) 



where g s is a universal function. From this it follows that 

/oo 
g s (E)dE 
-oo 

/oo 
g s (x)dx = CL- 1/v , (3.9) 
-oo 

where C is a constant independent of the system size. 
In order to perform the disorder averaging we consider 
40,000 samples for = 12, 10,000 for = 20, 2,000 for 
= 80, 1,000 for = 300, and 200 for = 1, 000. In 
each case the system consists of N z = 2 layers each with 
N,p states. We have also preliminary results for N z > 2. 
Since we exclusively consider systems with L x = L y we 
shall in the following use L — l^j2nN ( f > . 

Building on the work of Thouless and coworkersEao 
that relates g s (0) to o xx Andco has proposed that a 
similar relation should hold in a magnetic field, 



<J XX (E) = lim -—g s (E,L). 

L^-OO Z tl 



(3.10) 



This is however not true in any strict sense and one would 
only expect the left hand side of the above expression to 
be proportional to <j xx . However, in the absence of a 
magnetic field the root-mean-square leveL curvature can 
be related to the dissipative conduct anceO. If the scaling 
theorya of the quantum Hall effect is correct we would 
expect that o xx = (l/2)e 2 /h at the critical point and 
thus g s (0) = I /it. This seems to be consistent with 
what is found numerically for short-range scatterersLjEHI 
although thc^aage of the potential can change the value 
significantlyoHO A more rigojraj^panpmach would be 
to calculate the Chern numbersLSlIIIIla'tlfED, which con- 
firms the results of the scaling theory. In our calculations 
we find g s (0) ~ 0.2, (see Fig. ||) where calculations are 



shown with (thick solid lin e) o r without (thick dashed 
line) the sum over s in Eq. (3T). Clearly the results are 
markedly different. Also the density of stat es d iffer sig- 
nificantly without the infinite sum in Eq. (3.1) we find 
a de nsity of states that is no longer well described by 
Eq. (3.5). In Fig. |^ we also show results where the geo- 
metric mean exp(< ln|Ai?| > av ) has been used instead 
of |AS|. With the sum over s in Eq. (3T) the Thouless 
number is indicated as the thin solid line in Fig. ^ and 
without the sum as the thin dashed line. 



B. Results 

Before we turn to a discussion of the numerical results 
let us begin by looking at a few simple limits: 

Correlated disorder: In this case we can treat 
the N z — layer case straight forwardly. By Fourier- 
transforming along the z~ direction we remove the off- 
diagonal tunneling elements and obtain a matrix that is 
block diagonal with N$ x blocks, M„ . If we denote the 
wave- vector along the z-direction by — 2nn/L z , n — 
1, . . . N z , where N z = L z /a is the number of layers, we 
find that each of these blocks can be written: 



M„ = B — 2f cosAfl, 



(3.11) 



where B is the matrix describing the disorder in one layer. 
The presence of the 2tcosA:^I term will not affect the 
Thouless numbers since it is independent of the bound- 
ary conditions in the xy plane. Thus, except for some 
accidental degeneracies we should find that the Thouless 
numbers are a simple superposition of the the 1-layer re- 
sult displaced by 2t cos : 

g s {E) = J29 1 s - lweI (E-2tco8k* n ). (3.12) 

n 

Strong tunneling, uncorrelated disorder: Now we con- 
sider the case of uncorrelated disorder. In the limit 
where t tends to oo the tunneling completely dominates 
over the disorder and it is again advantageous to per- 
form a Fourier transform in the z— direction. In the limit 
t — > oo we can neglect the block-off-diagonal matrices and 
we again obtain a block-diagonal matrix with x 
blocks, M n . However, this time we find: 



M, 



— YB r 



2t cos 



(3.13) 



where B m is the matrix describing the disorder in the 
mth layer. Let us consider the two-layer case, N z — 
2. For p(E) we obtain two widely separated peaks each 
corresponding to a single layer with double the number 
of impurities at half the strength. Each of these peaks 
then have a width T = F/y/2. Hence, the integrated 
Thouless number, A(L), should increase by a factor of 
\J~2 compared to the single layer result. 
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Zero tunneling, uncorrelated disorder: If we set t = 
the Hamiltonian matrix again becomes block-diagonal 
and we should obtain results similar to the single layer 
case for large enough system sizes. The density of 
states, p(E), should remain unchanged and g s (E) ~ 
N z gl~ layel (E), since the total energy of states, N(E), 
is proportional to N z . 

We shall mainly be concerned with disorder that is not 
correlated between the layers (the top panel in Fig. [I]) 
but we shall briefly also discuss the case of correlated 
disorder (the bottom panel in Fig. The bulk of our 
results are shown in Fig. ^ where we display the den- 
sity of states, p(E), along with the Thouless number, g s , 
for a range of couplings, t = 0.0, 0.05, 0.15, 0.20, 0.25, 
between the two layers. In all cases the disorder was 
taken to be independent in the two layers. For the un- 
correlated disorder model we consider here, the effective 
correlation length of the disorder potential is the micro- 
scopic length I. We therefore expect that tunneling will 
be important when the localization length of decoupled 
layers exceeds £d r ~ £T/i, or, equivalently, £ exceeds l^ r . 
Accordingly the influence of tunneling on the density of 
states and especially on the Thouless numbers appears 
first near the center of the levels where the localization 
lengths are large. Well separated localized states appear 
for sufficiently large t/T. In the top row our results for 
two uncoupled layers, t = 0.0, are shown. We see that 
the Thouless number is exactly twice the value of a single 
layer (see Fig. |3|). For very weak tunneling the peak in 
g s becomes significantly broadened and for the sizes we 
have considered we do not observe two separate peaks un- 
til i ~ 0.1. As the tunneling between the two layers, t, is 
increased two peaks corresponding to the antisymmetric 
and symmetric state become apparent. For t = 0.25 we 
find that g s and the density of states essentially behave as 
the superposition of two independent peaks for the sym- 
metric and antisymmetric state in agreement with our 
considerations above. We also observe that the width 
of the density of states in the individual peaks approxi- 
mately obeys the relation T ~ ri_i ayer /\/2. Our results 
are in-good agreement with prior calculations by Ohtsuki 
et. alE3. In fig. ^ the position of the degenerate Landau 
levels in the absence of disorder is indicated as dashed 
lines in the panels for g s . For large t, level repulsion is 
clearly visible and the two extended state energies are 
further apart when disorder is included. 

Before discussing the scaling properties of the transi- 
tion we compare results for disorder that is independent 
in the two layers (uncorrelated disorder) or the same (cor- 
related disorder). In Figs. [|, || we show the density of 
states and g s for the case of uncorrelated disorder in two- 
layers coupled with a tunneling parameter of t = 0.1. The 
density of states shown in Fig. | is markedly broader than 
what we found for two uncoupled layers (top row Fig. |J). 
In Fig. ^| we show the Thouless numbers for t = 0.1, 
only for the largest sizes does it become clear that in the 
thermodynamic limit extended state energies exist at two 
discrete energies rather than across a band of finite width 



max - 0.20 for 
0.0) for two 



between low and high energy mobility edges. Note also 
that the peak value of the Thouless number is in this case 

= 0.1) ~ 0.22 whereas we found g 
the single layer case and <?™ ax ~ 0.35(< 
layers in the absence of any tunneling. 

The case of correlated disorder is very different. In 
Figs. [?], U we show our results for the density of states and 
g s , respectively. The disorder in this case is the same in 
the two layers and the tunneling parameter was, as above, 
taken to be t = 0.1. If we compare p(E) for uncorrelated 
and correlated disorder to the single layer results we see 
that p(E) is decreased at E = in both cases but more so 
for the case of uncorrelated disorder. This corresponds to 
a depletion of states at the center of the Landau band. 
The Thouless number, g s , is also significantly different 
for the case of correlated disorder. g™ ax is approximately 
0.28 for the smaller sizes before approaching a value of 
j s max ~ 0.20 for the largest systems. As expected the 
Thouless numbers are well described by a simple super 
position of single layer results Eq. ( 3.12 ). This is clearly 
not the case for the results in Fig. 

The difference between correlated and uncorrelated 
disorder is also reflected in the scaling of integrated Thou- 
less number, A(L) Eq. ( |3.9| ). In Fig. ^ we show results 
for the integrated Thouless number for thr ee different 
tunneling strengths. We see from Eq. (3.9) that A(L) 
according to the finite-size scaling form should behave as 
a power-law in L with exponent — 1/v. For the case of un- 
correlated disorder, the results in Fig. ^| yield 1/v ~ 0.43 
for t = 0.25 in agreement with previous results on sin- 
gle layer systems. Correlated disorder (plusses in Fig. ||) 
behaves in a similar way and we find 1/v ~ 0.44. How- 
ever, the case of uncorrelated disorder with weak tun- 
neling shows a marked crossover. For t — 0.1 we find 
1/v ~ 0.24. (Interestingly, numerical calculations for 
single-layer systems in the N = 1 Landau level show 
similar apparent enhancement of the localization length 
exponent.) For t — 0.05 we see that the smaller system 
sizes show the 1/v ~ 0.44 behavior expected for decou- 
pled systems before crossing over to a different power law. 
For small t and short-range potential correlations we do 
not expect tunneling to have any effect until the system 
size reaches ~ 1/t. If we fit only to points with L > 15 
we find, for t — 0.05, 1/v ~ 0.23 in very good agreement 
with the result for t = 0.1. For t = 0.15 and t = 0.2 wc 
find in both cases that the slope of A(L) increase with 
L, without saturating for the values of L available. This 
is consistent with the results shown in Fig. ^ where two 
separate peaks in g s clearly are visible at large L for these 
two values of t. 

We have also tried to analyze the Thouless numbers by 
integrating separately over the regions inside and outside 
the extended energies. This is difficult to do for small 
t since the extended energies cannot be located with a 
very high precision. Analyzing the integrated Thouless 
numbers separately for the two regions it is clear that 
the number of extended states between the two extended 
energies decreases significantly slower with L than for the 
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region outside the extended energies. 



IV. NETWORK MODEL 

We now proceed to discuss our results for a network 
model of the double-layer system. The network model 
was introduced by Chalker and CoddingtonO to take into 
account the corrections to percolative behavior that oc- 
cur when the correlation length diverges in the vicinity of 
the extended state energies in the middle of the Landau 
level. It is possible to map the network model for a sin- 
gle layer on to various spin modelsEjO and the calculated 
effective-correlation lengths can be used to estimate <j xx 
and a X yEl. The model that we use is essentially identical 
to one_±hat has been studied in previous work by Lee 
et. alE3. Each individual layer is represented by a sep- 
arate network model in the manner described in Ref. ^. 
A question now arises as to how to include tunneling 
between the layers. We follow Ref. [5?] and introduce a 
second saddle point along the straight paths in the orig- 
inal network model, coupling the two layers. A priori 
there are several ways to represent such a saddle point 
by a matrix. Since we want to model two physical layers 
and not pseudospins we make a slightly different choice 
than Ref. [57|. We take the inter-layer saddle point to be 
identical to the intra-layer saddle points, i.e. represented 
by the following matrix: 



_ /M° 
T '-{ Ml 



(4.1) 



with M™ given by: 

( e^" \ / cosh 







sinh 9t 
\ sinh 9 t cosh 9 t 





J<t>2 



(4.2) 



Here 9t is the parameter that controls the tunneling be- 
tween the two layers and n is the channel index counting 
the number of channels in each layer. We shall always 
take 9t to be constant. Random phases are included 
along the strai ght paths as in Ref. || described by the 
(j)'s in Eq. (4.2). The saddle points in the two layers are 
represented by identical matrices but with a different pa- 
rameter, 9, which we again take to be a constant and the 
same in the two layers. 



' cosh 9 a sinh 9 a 

cosh 9\, sinh&b 

sinh# a cosh# a 

sinh 9b cosh 9b , 



(4.3) 



Here the index a, b refers to the two layers. We shall 
always take 9 a — 9b since we are interested in modeling 
layers with equal density. 

Th e c hoice of the matrix coupling the two layers, 
Eq. (4.2), is by no means obvious. We could have 



used trigonometric functions instead of hyperbolics as in 
Ref. 57 thereby implying that in the picture where an in- 
dividual layer is represented by coupled lines of opposite 
going currents (see for instance Ref. |5^) the two layers are 
stacked in register with the currents going in the same di- 
rection in the two layers. One can also stack the layers out 
of register and the choice we have made corresponds in a 
certain sense to a mixture of these two choices. We be- 
lieve that none of these microscopic details should matter 
for the universal properties of the model and in particu- 
lar for the divergence of the correlation length at least for 
small 9 t . We have explicitly checked this for the results 
presented below in Fig. |l^, by repeating the calculation 
for several other choices of stackings and coupling matri- 
ces. For the small coupling of 9± = 0.05, used in Fig. [H], 
no dependence on the microscopic results was observed. 

We determine the correlation lengths associated with 
double-layer systems, described by the transfer matrices 
outlined above, by estimating the Lyaponov exponents. 
The positive Lyaponov exponents, \ l M {9), and their un- 
certainties are calculated following the method in Ref. |8l] 
for a range of values of 9 for fixed 9t- The correlation 
length is determined as the inverse of the smallest posi- 
tive Lyaponov exponent, X\ /I (9). 



Z M {6) = l/\ x M {9). 



(4.4) 



It is only necessary to calculate the positive Lyaponov 
exponents which saves considerable computing time. An 
additional check on the calculation can be done by cal- 
culating the first (and smallest in absolute value) of the 
negative Lyaponov exponents which should be the nega- 
tive of \\ji9jB. 

In general the correlation length, will be lim- 

ited by the width of the strip, M. The relevant quan- 
tity to study is therefore the reduced correlation length, 
Aa/(#) = £,m/M. An insulating region will be character- 
ized by Am(9) — > 0, a metallic one by Am(9) — > oo, or 
constant. Since Am(9) is dimensionless, standard finite- 
size scaling arguments predicts the scaling form: 



(4.5) 



where £ is the correlation length in the infinite system. 
As the critical energy, E C) is approached this correlation 
length diverges with the exponent v, 



£~\E-E C 



(4.6) 



The relation between 9 and the distance to—the critical 
energy, j c , can be determined approximatelyE3E.il for pos- 
itive 9, 



7 w lnsinhfl, 9 > 0, 



(4.7) 



implying that the critical 9 is given by sinh 9 C = 1, or 9 C = 
0.8814.... The relation Eq. (fO), allows us to rewrite 



the finite-size scaling relation, Eq. (4.5), in the following 
form, 
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Am(0) = giiM 1 '"), 



(4.8) 



which is the form we shall use in the analysis of the nu- 
merical results. 



A. Computational Method 

We perform the calculations in a cylinder geometry 
imposing periodic boundary conditions in the transverse 
direction. The system is made invariant under a rota- 
tion by 90° by alternating transfer matrices with 9 and 
sinh 9' — 1/ sinh 9 as described in Ref . ^. We denote the 
width of each of the two layers by M and consider sys- 
tems with M ranging from 4 to 128. For M = 4, 8 we 
generate 2 x 10 6 transfer matrices and for the remaining 
widths 2 x 10 5 . We obtain sets of data by fixing 9 t and 
approaching the critical point by varying 9 and thereby 7. 
As a check on our calculations we set 9t to zero and were 
able to reproduce the single-layer results from Ref. || to 
within statistical errors. The integral factor sometimes 
introduced in the definition of the number of layers was 
chosen so that this would be the case. 



B. Results 

Our main results on the network model for the double- 
layer system is presented in Fig. |o|. For lattice sizes 
ranging from M = 4 to M — 128 we have calculated 
the reduced correlation length, Am (9) . Since we want to 
view the tunneling between the two layers, described by 
Of, as a small perturbation we take this parameter to be 
very small and constant, 9 t =0.05. Roughly we have the 
relations tanh# t = t, where t is the parameter describ- 
ing the tun neling in the exact diagonalization studies in 
Section III. Note, that we get t and not 2t since we do 



not have periodic boundary conditions between the two 
layers as we had in the Section 111. We then vary the 
intra-layer coupling, 9, and plot the results as a function 
of 7. As clearly seen in Fg. 10| we do not reach a scaling 
regime until the width of the strips, M , exceeds 1/t ~ 20, 
as expected from the discussion in Section H We expect 
the extended state energi(es) to be located at 7 ~ since 
we have taken the intra-layer coupling, 9, to be the same 
in the two layers. For widths larger than M=20 we ob- 
serve that the reduced correlation length becomes inde- 
pendent of the width, M, at 7 = as expected. For the 
sizes considered we do not see any signs of a splitting of 
the two extended state energies which should have been 
of the order of A7 ~ 2t ~ 29 1 , based on the simple tight- 
binding picture, and therefore clearly visible in Fig. |lC| . 

Given these observations we therefore perform a scal- 
ing analysis under the assumption that both of the two 
expected extended state energies are to be found at the 
same critical energy 7 C = 0. Since we do not expect scal- 
ing to be obeyed for M -C 20 we only include widths 
M ^> 20. Testing the scaling analysis is now a simple 



matter of rescaling the x-axis in Fig. [To] by an amount 
M x / V for the different widths. The result is shown in 
Fig. |l[ Clearly very good scaling is found for the chosen 
value of the correlation length exponent v — 14/3. Since 
we do not have a large number of numerical data avail- 
able to determine the exponent, we can only test if a 
given value of v gives good scaling. We tried v = 7/3 and 
v = 11/3 in both cases we found scaling that visibly was 
much worse than what is shown in Fig. pd| with v = 14/3. 
We therefore conclude that the apparent doubling of the 
correlation length exponent is not in disagreement with 
the numerical results obtained from the simplified net- 
work model for the double-layer system. 

We now wish to make a few comments on the applica- 
bility of the network model to the real physical system. 
We take the view that the starting point for the network 
models is a percolation path close to a critical energy. 
Quantum tunneling at the saddle points then intro- 
duced as a small correction to the physicsE£l. We believe 
that the double-layer network model only describes the 
physics of the coupled quantum wells in detail, in the 
limit where the inter-layer tunneling parameter tends to 
zero, 9t — > 0. Our argument goes as follows: Imagine we 
wanted to build a network model describing a strongly 



coupled double-layer system. As we saw in Section III 
the symmetric and antisymmetric states are then widely 
split, by an amount As^s ~ 4i. As t becomes very large 
we should therefore model the system as two weakly cou- 
pled single layer networks, describing the symmetric and 
antisymmetric states. These two networks should then 
be modeled using different interlayer couplings 9\ and 
#2 corresponding to the symmetric and antisymmetric 
states. In the limit t — ► 00 the two networks become com- 
pletely decoupled (and \6i — 9 2 \ — ► 00) and we should thus 
set 9 t — 0. Related ideas were proposed for a network 
model in Ref. |35| in order to describe Landau level mix- 
ing and related ideas have also previously been discussed 
in Ref. [58]. In our model, as we have described it above, 
the two networks correspond to the two layers in physical 
space. As the intra-layer coupling is increased we do not 
form symmetric and antisymmetric states, but rather the 
paths become localized in orbits between the two layers. 
This is clearly seen in Fig. 12 where we plot results for 
9 t = 0.5 and 9 t — 2.0. Clearly, the reduced correlation 
length is growing slower than M, and for 9 t = 2.0 it ap- 
pears that it is independent of M and also of 7, consistent 
with our expectation that the orbits all are localized be- 
tween the two layers on a length scale of order 1. This 
is supported by the observation that for Q t = 2.0, Am (0) 
decreases with M roughly as M _1 indicating a constant 
correlation length. For intermediate couplings it is possi- 
ble that one could see structure resembling the symmetric 
and antisymmetric state but we believe that the double- 
layer network models only describe the correct physics 
of the real double-layer systems in the limit where the 
inter-layer tunneling parameter, 9 tl is but a small pertur- 
bation. Hence, we believe that the double-layer network 
model do not describe the microscopic physics in detail 
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for intermediate tunneling parameters. 

V. SUMMARY AND DISCUSSION 

We have shown that in the strong magnetic field limit 
non-interacting double-layer electron systems in which 
interlayer tunneling occurs have two extended state ener- 
gies for each orbital Landau level. Numerical results show 
that, for large tunneling amplitude, Landau levels asso- 
ciated with subbands which are symmetric and antisym- 
metric combinations of isolated layer states are weakly 
mixed by disorder. Localization properties within the 
Landau levels of the symmetric and antisymmetric sub- 
bands are similar to those for a single 2D electron layer 
and in particular appear to have a correlation length ex- 
ponent, v, identical to what was found for an isolated 
layer to within numerical precision. For smaller values of 
the tunneling amplitude where symmetric and antisym- 
metric subbands are not well developed in the density 
of states, numerical results still appear to show that ex- 
tended states occur at only two energies which are split 
by an amount somewhat larger than the splitting of sym- 
metric and antisymmetric Landau levels in the absence 
of disorder. We cannot exclude the possibility that a fi- 
nite amount of tunneling is necessary to split the two 
extended state energies, although the weight of available 
evidence appears to suggest the contrary. However, nu- 
merical values of the localization length are much larger 
than in the limits of either strictly zero tunneling or large 
tunneling. This is especially true in the energy interval 
between the two extended state energies. Over the range 
of system sizes accessible for numerical studies the expo- 
nent for the diverging localization length appears to be 
approximately twice as large as for the case of isolated 
layers. We have compared these microscopic calculations 
with network models of double-layer systems. The net- 
work model and microscopic numerical results differ qual- 
itatively. In the network model case there is no evidence 
for two discrete critical energies at which extended states 
occur. We conclude from our study that, in contrast with 
the single layer case, network models do not generically 
give reliable results for the strong magnetic field localiza- 
tion properties of double-layer systems. 

We believe that our results have important impli- 
cations for the integer quantum Hall effect in three- 
dimensional electron systems. We comment here only 
on the case where the band width along the field di- 
rection is smaller than the Landau level separation and 
the quantum Hall effect has the best chance of occur- 
ring. It is important to realize that the physics of this 
extreme strong field regime is qualitatively different from 
the more usual three-dimenakinal case where many dif- 
ferent Landau 'tubes' crossEEl the Fermi energy. (The 
physical systems we have in mind are multiple quantum 
well (MQW) systems, like those studied experimentally 
by Stormer et al., with weak barriers between the wells.) 



In particular, just as in the single layer integer quantum 
Hall case, we can argue that disorder can never result in 
the localization of all states. This point is perhaps made 
most elegantly using the-jtopological picture of the inte- 
ger quantum Hall cffcctED. In the absence of disorder an 
elementary calculation shows that the Hall conductance 
in e 2 /h units, and hence the sum of the Chern numbers of 
all states, is equal to the number of layers in the MQW. 
As the states evolve adiabatically with disorder, the sum 
of all the Chern numbers o£-all states in the (energy range 
of interest) cannot changeO. Since only extended states 
can have non-zero Chern numbers, it is impossible to lo- 
calize all states. The situation is closely analogous to the 
quantum Hall effect in a single two-dimensional electron 
gas in a magnetic field when Landau level separations 
become small, since each Landau level contributes I to 
the Chern number sum and localization is possible only 
by mixing states with different Chern number which are 
at energies well away from the bottom of the band of the 
host semiconductor. 

In the absence of disorder the states in the energy range 
of interest in the MQW consist of a set of macroscopi- 
cally degenerate Landau levels, labeled by wave-vectors 
and split by an amount proportional to the interlayer 
hopping amplitude. For small T/t, states with a given 
wave- vector will be weakly coupled and a single extended 
state energy will exist for each Landau level. As the dis- 
order strength increases, our numerical results for the two 
layer case suggest that the extended state energies will re- 
main separate but localization lengths will increase sub- 
stantially, except at energies above the highest extended 
state energy and below the lowest extended state energy. 
(However, we know of no general argument which forbids 
either the collapse of the extended states toward a single 
energy or, in the other extreme, the development of a 
band of energies over which states are extended.) In the 
limit of an infinite number of layers, the energy separa- 
tion between extended states will approach zero but the 
system will still, strictly speaking, not be metallic siuee 
almost all states in any range will still be localizeoS. 
Our numerical results suggest the possibility that critical 
exponents for localization lengths diverging between in- 
termediate energy extended state could be different from 
the critical exponents for single-layer systems. 

It is interesting to consider whether or not the metallic 
phase of three-dimensional systems in a strong magnetic 
field suggested by the work in Refs. ||^J5|,|5^, can be rec- 
onciled with the expectation from integer quantum Hall 
theory and from the present calculations of a discrete set 
of extended state energies for any finite number of lay- 
ers. In order for these two pictures to be compatible the 
localization lengths at energies between the lowest and 
highest extended state energies would have to increase 
with the number of layers and diverge in the limit of infi- 
nite layer numbers or equivalently the limit of small sep- 
arations between extended state energies. Incidentally, 
if we assume that each additional quantum well leads 
to another extended state energy without changing the 







associated correlation length exponents as this energy is 
approached, it is not clear how to explain the large differ- 
ence between the exponent v found at the mobility edge 
in Refs. |5l],|52|,^5] and the two-dimensional exponent. We 
believe that the tendency toward an apparent metallic 
phase in multiple quantum well systems in the limit of 
large layer numbers or small interlayer hopping is ex- 
tremely closely connected with the disappearance of the 
quantum Hall effect in a high-mobility two-dimensional 
electron system in the limit of weak magnetic fields since 
disorder in both cases permits only mixing of Landau 
bands carrying the same unit total Chern number. Ex- 
isting experiments on the integer quantum Hall effect 
in MQW systems have observed a quantized Hall effect 
only at Fermi energies above the highest energy extended 
state where electrons are well localized. We hope that 
the present paper will motivate new attempts to study 
the physics of the quantum Hall effect in double quan- 
tum well and MQW systems at Fermi energies between 
extended state energies. 
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FIG. 3. The Thouless number, j s , as a function of energy 
for a single layer system with N$ = 72, and n p = 5N$. Cal- 
culations are shown using the elliptic theta functions (thick 
solid line) as well as without them (thick dashed line). We 
also show results for g s calculated using a geometric mean and 
the theta-functions (thin solid line) and without them (thin 
dashed line). 
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FIG. 1. The two different forms of disorder. 




FIG. 2. The flow diagram in the presence of tunneling, 
from Ref. §t|. 
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FIG. 4. The density of states, p(E), and the Thouless num- 
ber, g s , for a range of couplings t between the two layers. 
n p — 3N,/, random pinning centers, modeled as random sign 
(5-functions, were used in the calculation. The position of the 
pinning centers was taken to be random and independent in 
the two layers, corresponding to uncorrelated disorder. In all 
figures the curves corresponds to iV^ = 12, 20, 80, 300, 1000. 
The energy E is in the units of Tscba- The dashed lines 
indicates the position of the extended state energies in the 
absence of disorder. 
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FIG. 5. The density of states for a coupling t = 0.1 between 
the two layers. n v — 3N,/, random pinning centers, modeled 
as random sign ^-functions, were used in the calculation. The 
position of the pinning centers was taken to be random and 
independent in the two layers. The two dashed lines indi- 
cates the position of the two degenerate Landau levels in the 
absence of any disorder. 



FIG. 7. The density of states for a coupling t = 0.1 between 
the two layers. n v — 2>N$ random pinning centers, modeled 
as random sign 5-functions, were used in the calculation. The 
position of the pinning centers was taken to be random but the 
same in the two layers, corresponding to a form of correlated 
disorder. The two dashed lines indicates the position of the 
two degenerate Landau levels in the absence of any disorder. 
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FIG. 6. The Thouless number, g 3 , for a coupling t = 0.1 
between the two layers. n p = 3iV^, random pinning centers, 
modeled as random sign 5-functions, were used in the calcu- 
lation. The position of the pinning centers was taken to be 
random and independent in the two layers. The two dashed 
lines indicates the position of the two states in the absence of 
any disorder. The results corresponds to the density of states 
in Fig. | 



FIG. 8. The Thouless number, g s , for a coupling t = 0.1 
between the two layers. n p = 3A r </ ) random pinning centers, 
modeled as random sign 5-functions, were used in the calcu- 
lation. The position of the pinning centers was taken to be 
random and but the same in the two layers, corresponding to 
a form of correlated disorder. The two dashed lines indicates 
the position of the two degenerate Landau levels in the ab- 
sence of any disorder. The results corresponds to the density 
of states in Fig. |^ 
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FIG. 9. The integrated Thouless number as a function of FIG. 11. Scaling collapse of the reduced correlation length 

L, for tunneling strenghts of t — 0.05, 0.1, 0.25. The plusses for the double layer network model with a tunneling param- 
denote results for t = 0.1 with correlated disorder. eter of 6 t = 0.05. 




FIG. 10. The reduced correlation length as a function of FIG. 12. The reduced correlation length for two different 

the energy 7 = lnsinhf? for the double layer network model values of t = 0.5 and 2.0. the results are for for the double 
with a tunneling parameter of t = 0.05. layer network model. 
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